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NEUTRAL  ENTRAINMENT  DEMONSTRATION  IN  A  XENON  FRC  THRUSTER  EXPERIMENT 

David  Kirtley,  George  Votroubek,  Anthony  Pancotti,  Jim  Pihl 
MSNWLLC,  Redmond,  WA. 

John  Slough 

The  University  of  Washington,  Seattle,  WA 

ABSTRACT 

The  primary  concern  with  operating  a  fully-ionized,  plasma-based  electric  propulsion  system  on 
lightweight  gasses  is  that  the  per-ion  plasma  formation  and  ejection  losses  are  typically  much  larger  than 
the  energy  of  directed  lightweight  propellant.  A  new  innovation,  the  PROTEAN,  adds  a  secondary  Neutral 
Entrainment  (NE)  stage  that,  if  succesfull,  can  reduce  the  plasma  formation  losses.  This  NE  stage  injects 
neutral  propellant  into  the  plasma  path  after  ejection  from  an  electromagnetic  Field  Reversed 
Configuration  (FRC)  thruster,  such  as  the  1  kW  Electromagnetic  Plasmoid  Thruster  (EMPT).  The  FRC 
ingests  these  neutral  particles,  accelerating  them  up  to  high  velocities  through  resonant  charge  exchange 
collisions.  The  NE  stage  then  highly  efficiently  adds  further  kinetic  energy  to  the  engorged  plasmoid 
through  Peristaltic  Dynamic  Acceleration,  in  which  sequenced  theta-pinch  coils  are  actuated  to  provide  a 
large  magnetic  pressure  gradient.  In  total,  this  system  should  allow  for  the  addition  of  active  propellant 
mass  without  significant  further  ionization  or  plasma  frozen  flow  losses.  This  paper  presents  the 
experimental  and  MHD  modeling  results  of  the  Phase  I  PROTEAN  program  to  demonstrate  Neutral 
Entrainment  in  a  5  Joule  conical  experiment.  Radial  and  temporal  plasma  profiles  were  measured  using 
external  and  internal  magnetic  field  probes  and  downstream  time-of-flight  langmuir  probes  for  the 
formation,  acceleration,  and  entrainment  of  a  xenon  FRC. 

INTRODUCTION 

Electric  Propulsion  (EP)  systems  accelerate  charged  particles  by  the  application  of  electrostatic  and 
electrodynamic  (magnetic)  forces.  In  order  to  accelerate  a  propellant  it  must  be  charged  and  in  typical 
electric  propulsion  systems  it  must  undergo  an 
ionizing  collision  event.  The  primary  loss 
mechanism  for  all  fundamental  non-thermal 
electric  propulsion  devices  and  technologies  is 
plasma  formation  loss  [1],  In  a  realistic  EP  system, 
it  takes  100-500  eV/per  atom  to  take  one  neutral 
propellant  particle  and  eject  one  ionized  particle 
not  counting  kinetic  energy  addition.  All  plasma- 
based  electrostatic  or  electromagnetic  propulsion 
systems  require  an  ionized  particle  to  be  able  to 
accelerate  propellant.  Fundamentally,  the  energy 
cost  of  an  ionizing  collision  is  quite  small,  on  the 
order  of  10-20  eV.  This  would  roughly  correspond 
with  an  accelerating  voltage  of  10-20V  for  50% 
efficiency.  However,  this  is  not  a  realistic 
assessment  of  the  total  energy  required  for 
ionization.  In  real  plasma,  energy  losses  also 
include  recombination,  excitation  (radiation),  and 
polarization.  These  effects  will  tend  to  dominate  the  effect  of  the  ionizing  collision,  specifically  at  low 
energies.  Figure  1  shows  a  compilation  of  energy  loss  mechanisms  for  an  ideal  plasma  in  argon  and 
oxygen  which  have  minimum  energies  of  ~14  and  ~16  eV/ion.  It’s  clear  that  for  electron  temperatures 
below  50  eV,  the  minimum  total  ionization  potential  is  much  greater  than  14-16  V.  Neutral  Entrainment 
attempts  to  minimize  the  largest  energy  loss  of  electromagnetic  thrusters  by  effectively  accelerating 
neutral  propellant,  without  net  ionization. 


Te(V) 

Figure  1.  Ionization  energy  cost,  including 
excitation,  as  a  function  of  electron 
temperature  [2]. 
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BACKGROUND  OF  PLASMA  SOURCE 


It  is  believed  that  a  neutral  entrainment  system 
can  operated  with  a  variety  of  plama  sources,  so  long  as 
that  plasma  is  collisional  and  electomagnetically 
accelerated.  For  this  experiment,  a  unique  plasma 
source,  a  Field  Reverse  Configuration  (FRC),  was  used 
which  has  closed  field  magnetic  topology  and  relatively 
high  plasma  density.  A  FRC  plasma  consists  of  a  closed 
field  line,  fully  magnetized  plasma  with  a  large  azimthal 
current  (as  in  Figure  2).  The  plasma  diamagnetic  current  Figure  2.  Cross  section  of  a  cylindrical 

is  generated  opposite  to  the  applied  external  field,  field  reversed  configuration  plasma, 

creating  a  high-pressure  configuration  [3],  Typical  FRC 
plasmas  are  formed  with  fast  (few  ps),  large  (10’s  of  kA) 
inductive  discharges. 

The  generation  of  FRCs  using  rotating  magnetic 
fields  (RMF)  is  well  established  [4],  The  transverse  rotating 
magnetic  field  is  produced  by  external  antennas  as  indicated 
in  Figure  3.  The  first  pair  of  loop  antennas,  carrying  a 
sinusoidally  occillating  current,  produce  a  uniform  magnetic 
field  transverse  to  the  coil  axis  as  shown  in  black.  A  second 
loop  pair  produce  a  transverse  magnetic  field  perpendicular 
to  the  first  pair  (blue)  and  are  driven  90°  out  of  phase  with 
respect  to  the  first  pair.  The  composite  magnetic  field  then 
forms  a  rotating,  transverse  magnetic  field  of  constant 
magnitude.  An  electron  magnetized  to  the  field  lines  will  rotate 
azimuthally  on  average  at  the  same  frequency  as  the 
antennas.  An  axial  bias  field  is  applied  and,  if  the  magnetude 
of  the  RMF  induced  Je  is  suffient  to  reverse  this  field,  a  true 
field  reversed  plasma  is  formed  as  in  Figure  4.  For  an  axially- 
distributed  RMF  system  it  is  found  that  the  azimuthal  current 
is  quite  simply: 

Jg  =  en/or 

The  PROTEAN  source  forms  a  RMF  in  a  conical, 
gradient  magnetic  bias  field.  The  large  Je  is  driven  in  a  conical 
field  with  a  radial  magnetic  field.  The  fully-formed  FRC  is  then 
accelerated  axially  by  the  resultant  JexBr  force.  The  azimuthal 
current  continues  to  be  generated  as  the  FRC  moves 
downstream  and  the  FRC  accelerates  throughout  the  entire 
cone.  Finally,  as  the  FRC  expands  through  the  conical  section 
and  beyond  the  exit  of  the  cone,  the  thermal  energy  of  the 
FRC  is  converted  into  axial  velocity  as  in  a  typical  magnetic  expansion,  resulting  in  a  low  temperature, 
high  velocity,  and  compact  plume.  The  key  parameter  that  will  be  discussed  is  that  the  FRC  current 
generates  a  self-magnetic  field  which  supports  a  high  internal  plasma  pressure.  This  pressure  is  in 
equilibrium  with  the  external  field  and  allows  for  magnetic  compression  and  acceleration,  as  will  now  be 
described. 


O  O 


Figure  3.  Schematic  of  an  RMF. 


Figure  4.  Schematic  of  a  Rotating 
Magnetic  Field  creating  an  FRC. 
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THE  PHYSICS  OF  NEUTRAL  ENTRAINMENT 


The  fundamental  concept  in  Neutral 
Entrainment  (NE)  is  that  a  FRC  will  ingest  large 
quantities  of  neutral  gas  through  charge  exchange 


0  5  10  15  20  25  30 

Temperature  [eY] 

collisions,  not  ionization.  If  one  accelerates  an  FRC 
while  providing  upstream  neutral  gas,  velocity  can  be 
specified  and  propellant  mass  added  with  very  high 
efficiency.  In  recent  FRC  experiments  it  was  found 
that  accelerating  an  FRC  into  and  through  a  large 
neutral  population  was  beneficial  to  the  FRC  [5,  6].  It  t 


Figure  5.  Neutral  Entrainment  ELF  Thruster. 

i  gained  mass  and  became  more  stable. 


Neutral  entrainment  is  believed  to  occur  as  follows.  A  collisional,  ionized  propellant  is  accelerated 
in  a  gradient  magnetic  field.  In  this  specific  case  it  will  be  in  the  form  of  a  closed  field,  isolated  FRC  that  is 
translating  in  a  constant  bias  magnetic  field.  A  slow,  cold  neutral  gas  population  is  introduced  in  front  of 
the  translating  plasma.  As  the  plasma  ions  collide  with  the  neutral  particles  a  positive  ion  will  interact  with 
an  atom  so  as  to  capture  a  valence  electron,  resulting  in  a  transfer  of  the  electron  from  the  atom  to  the 
ion.  If  this  occurs  between  populations  with  the  same  mass,  then  the  charge  transfer  occurs  with  no 
change  in  radial  or  axial  momentum  and  no  decrease  in  kinetic  energy. 

After  this  collision  the  resulting  fast  moving  neutral  particle  proceeds  out  of  the  thruster 
unhindered  by  magnetic  fields  with  its  original  momentum.  Meanwhile,  the  slow-moving  charged  particle 
is  still  confined  in  an  accelerating  magnetic  field  and  is  subsequently  accelerated  to  high  velocity  and 
applying  a  thrust  force  to  the  spacecraft  (via  magnetic  field  coils).  In  total,  there  are  two  fast  moving, 
propellant  particles  but  only  one  ionization  loss,  effectively  reducing  the  total  ionization  loss  per  particle  in 
half.  Depending  on  the  collision  frequency,  velocity,  and  density  of  the  plasma  and  neutral  populations 
this  process  will  be  repeated  throughout  an  FRC  and  possibly  several  times.  It  is  expected  that  there  is  a 
maximum  of  10  neutral  entrainments  per  particle  possible  with  a  typical  NE  geometry  (simply  due  to  the 
collision  frequency  for  a  typical  geometry)  reducing  an  ionization  cost  of  40  eV/ ion  to  as  low  as  5  eV/ion, 
well  below  even  the  theoretical  minima.  Shown  in  Figure  5  is  an  MHD  simulation  of  an  accelerated  FRC. 
The  key  parameter  in  this  analysis  is  the  ratio  of  ionization  to  charge  exchange  collision  cross  section.  As 
is  well  known,  the  charge  exchange  collisions  are  much  more  likely  than  the  ionization  collisions  [2], 


Shown  in  Figure  6  are  measured  charge  exchange  and  ionization  cross  sections  for  xenon.  For 
the  energies  of  interest  charge  exchange  is  30+  times  as  likely  as  ionization.  Finally,  even  if  there  is  a 
10%  population  of  ionized  propellant  particles,  they  are  accelerated  with  the  propellant  and  are  not 
considered  to  be  a  loss.  Key  design  parameters  are  the  mean  free  path  for  collisions  for  both  charge 
exchange  and  neutral-neutral  collisions  for  the  length  of  the  neutral  entrainment  region  which  determine 
the  maximum  number  of  entrainments.  Additionally,  a  unique  benefit  of  the  plasma  is  that  they  typically 
do  not  relax  to  Maxwellian  energy  distributions  and  so  the  ionization  potential  of  the  high-energy  electron 
tails  is  not  realized  [8]. 
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The  Plasma  Science  and  Innovation 
Center  at  the  University  of  Washington  has 
developed  a  two-fuid  model  of  interaction  plasma 
and  neutral  fluids  [9].  This  model  is  based  on  the 
SEL-HiFi  code  [10],  a  3D  MHD  code  with 
advanced  features  that  are  ideal  for  modeling  this 
complex  interaction.  To  this  code  the  capability  to 
investigate  a  neutral  fluid  with  resonant  charge 
exchange,  electron-impact  ionization,  and 
radiative  recombination  reactions  was  added. 
Figure  7  show  results  from  a  propagating  FRC 
interacting  with  a  background  neutral  gas.  In 
these  cases,  the  total  plasma  mass  of  the  FRC 
and  neutral  mass  of  the  interaction  region  are 
equal.  Additionally,  the  neutral  mass  is  centered 
on  axis  as  in  the  experiment  and  extends  for  50 
cm  downstream.  In  these  simulations  a  30 
microgram  FRC  is  colliding  with  a  30  microgram 
neon  region  after  being  initialized  to  26  km/s  axial 
velocity.  Figure  7  shows  the  plasma  pressure 
increase  as  the  FRC  interactions  with  the  neutral 
pressure.  Finally,  as  one  would  expect,  there  is  an  increase  in  diamagnetism  and  a  slowing  of  the  FRC  as 
it  passes  through  the  gas. 

Figure  8  shows  the  temporal  profile  of  the  total  momentum  components  of  the  FRC  as  well  as  the 
neutral  population.  Again,  this  is  for  a  case  in  which  the  neutral  mass  was  equal  to  the  plasma  mass. 
Quite  obvious  from  these  results  is  that  both  charge  exchange  and  ionization  play  a  critical  role  in  these 
processes  and  that  momentum  can  successfully  be  transferred  to  the  neutral  gas  as  desired.  For  an 
equal  plasma  to  neutral  mass,  approximately  50%  of  the  momentum  of  the  plasma  is  transferred  to  the 
plasma,  with  a  majority  of  that  momentum  transfer  coming  from  charge  exchange.  In  these  cases,  it 
appears  that  there  is  significant  ionization  of  the  background  plasma  although  that  energy  appears  to 
come  from  the  initial  thermal  energy  of  the  plasma  rather  than  the  kinetic  energy,  exactly  as  desired. 
While  this  not  a  loss  mechanism  in  neutral  entrainment,  as  the  new  plasma  will  be  accelerated  with  the 
initial  FRC  plasma,  it  is  less  neutrals  entrained  overall,  lessening  the  efficiency  and  T/P  gain.  It  was  found 
that  this  ionization  ratio  is  a  direct  consequence  of  the  plasma  temperature  at  the  exit  of  the  thruster.  For 
5  eV  total  temperatures  (similar  to  experiments)  and  high  downstream  neutral  density,  it  was  found  that 
as  low  as  20%  of  neutrals  are  ionized  while  80%  of  the  plasma  momentum  is  transferred  to  the  neutrals. 


Figure  6.  Measured  xenon  Election  Impact 
Ionization,  Electron  Ionization  Excitation,  and 
Charge  Exchange  cross  section  [6,  7]. 

t=9  (as  t=17  jas 


Figure  7.  Pressure  psuedocolor  plots  for 
baseline  ELF  simulation.  Neutral  entrainment  is 
shown  by  white  arrows,  contour  is  normalized 
neutral  pressure. 


Figure  8.  Temporal  profiles  of  momentum-related  quantities  for  ELF  baseline  simulation 
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BACKGROUND-  PERISTALTIC  DYNAMIC  ACCELERATION  (PDA) 

Peristaltic  dynamic  acceleration  is  a  concept  that  has  been  demonstrated  in  several  FRC 
programs  at  MSNW  LLC  [11].  An  FRC  which  is  confined  within  a  cylindrical  magnetic  field  coil  can  be 
compressed  or  accelerated  by  that  axial  magnetic  field.  The  key  relation  to  dynamic  magnetic 
acceleration  is  Equation  1,  which  relates  the  internal  plasmoid  density  and  temperature  to  the  external 
magnetic  pressure.  This  is  for  a  magnetized,  fully  ionized  plasma  in  radial  equilibrium  with  the  external 
field.  The  magnetic  Beta  shows  the  amount  of  compression  of  the  FRC  relative  to  the  vacuum  condition, 
where  xs  is  the  normalized  magnetic  separatix  radius. 


P0  =  n0kT  = 
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Finally,  for  a  plasma  in  radial  equilibrium  a  relationship  for  magnetic  flux  conservation  (Equation 
2)  equates  the  magnet  flux  during  a  compression  cycle  to  flux  after  the  plasmoid  has  been  ejected.  This 
can  be  directly  translated  into  velocity  and  is  demonstrated  in  the  laboratory  at  very  high  efficiencies. 


n  n  1 

(Ar  -  An )L„  — ^  A„L„  =-Mvl 

2fto  P  P  2ft,  p  p  2  p 


(4) 


The  acceleration  sequence  can  then  be  represented  by  Figure  5  which  shows  the  magnetic  flux 
contours  for  a  translating  FRC  in  a  gradient  field.  In  this  simulation  the  FRC  is  formed  earlier  and  then 
translated  in  the  acceleration  section.  Three  pulsed  magnet  field  coils  are  triggered  in  sequence  (typical 
accuracies  <100  ns  required)  as  the  FRC  translates.  This  creates  a  net  magnetic  field  gradient  behind  the 
plasmoid  accelerating  it  to  high  speeds.  Additionally,  a  background  field  can  be  applied  to  act  as  a  guide 
field,  magnetic  wall  insulation,  and  to  prevent  reverse  voltages  from  being  applied  to  the  charging 
capacitors.  At  this  point  the  IGBT  switches  can  be  opened,  recovering  any  stored  inductive  energy  back 
to  the  capacitors.  They  are  then  recharged  for  the  next  discharge.  The  Inductive  Plasmoid  Accelerator 
used  sequenced  magnetic  field  coils  to  accelerate  a  theta-pinch  formed  FRC  to  greater  than  300  km/s 
and  required  timing  precision  of  less  than  100  ns. 

The  MHD  code  MOQUI  was  used  to  model  the  dynamic  acceleration  of  a  set  of  post  ejection 
acceleration  coils.  In  MOQUI,  an  FRC  was  formed  of  similar  densities  and  temperatures  to  the  RMF- 
formed  FRCs  in  EMPT.  It  was  then  allowed  to  travel  into  a  conical  section  with  a  single  active  magnetic 
field  coil.  In  Figure  9,  the  FRC  was  allowed  to  translate  freely  within  the  flux  conserving  region  and 
appears  to  be  a  translating,  well-formed  FRC.  In  Figure  10,  the  downstream  acceleration  coil  was 
activated,  compressing  and  accelerating  the  plasmoid.  From  the  figures  it  is  clearly  evident  that  the 
acceleration  process  occurs  as  expected,  with  the  plasmoid  increasing  in  velocity  and  the  primary  mass 
compressing  radially  as  it  accelerates  axially. 
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Figure  9.  An  FRC  translating  without  dynamic  Figure  10.  An  FRC  translating  with  dynamic 
acceleration.  acceleration. 
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EXPERIMENTAL  DESIGN 


The  experiment  design  is  optimized  for  a  30  kW  architecture  and  utilizes  a  conical  quartz 
insulator,  28  cm  long  and  14  cm  in  outer,  large  diameter.  The  experiment  has  a  total  of  ten  accelerator 
coils  axially  along  the  entire  length  of  the  experiment  setup  as  shown  in  Figure  12.  These  magnetic  field 
coils  integrate  both  the  bias  and  accelerator  fields  for  the  RMF  as  well  as  the  accelerator  section.  In  this 
way  the  optimal  balance  of  FRC  formation,  acceleration,  and  entrainment  can  be  determined  empirically 
without  a  preconceived  reliance  on  the  early  model  as  timing  can  easily  be  prescribed.  The  insulator  is  a 
12°  quartz  cone.  There  are  five  sets  of  external  magnetic  field  and  magnetic  flux  diagnostics  along  the 
length  of  the  experiment.  In  this  design  there  are  no  external  flux  conservers  as  there  have  been  in 
previous  thruster  manifestations.  In  the  PROTEAN  the  flux  conserving  boundary  is  established  with  the 
aforementioned  multi-turn  field  coils.  These  twenty  turn  coils  have  an  inductance  ranging  from  3 
microhenry  to  30  microhenry.  The  RMF  antennas  are  scaled  to  5  kW  and  placed  externally  to  the 
accelerator  coils.  One  of  the  critical  innovations  to  this  concept  is  the  idea  that  FRC  formation  and 
ejection  no  longer  has  to  be  a  passive  process.  The  FRC  plasma  source  plasma  can  be  injected  into  the 
entrainment  and  acceleration  section  at  the  optimal  density  and  plasma  temperature,  not  necessarily  the 
maximum  acceleration  pressure.  In  this  way  the  optimal  plasma  conditions  for  neutral  entrainment  can  be 
prescribed. 


GAS  INJECTOR 

One  of  the  most  complicated  aspect  of  the  PROTEAN 
experimental  setup  is  the  steady  injection  of  neutral  propellant 
to  be  entrained.  In  previous  Neutral  Entrainment  research  it 
was  clearly  shown  that  the  geometric  and  temporal  distribution 
of  neutrals  is  critical  to  successful  entrainment  [12].  The 
neutrals  must  be  located  near  the  axis  to  prevent  interaction 
with  the  quartz  walls  which  leads  to  compression,  and  thus 
ionization  of  the  neutrals.  Additionally,  the  neutrals  must  be 
placed  downstream  of  the  FRC  formation  section,  as  the  RMF 
is  highly  efficient  at  ionizing  neutral  gas,  resulting  in  complete 
40  eV/ion  ionization  of  the  propellant.  Therefore,  the  special 
neutral  injector  was  designed  as  shown  in  Figure  11  and 
Figure  13.  In  this  design  there  are  two  propellant  injection 
manifolds  that  can  provide  steady  gas  injection  independently.  Figure  11.  Photograph  of  the  high 
The  inner  manifold  is  a  Supersonic  Molecular  Beam  Injector  speed  SMBI. 

(SMBI)  that  uses  a  400:1  expansion  nozzle  to  accelerate  and  cool  a  neutral  gas  to  extremely  high 
velocities.  This  neutral  gas  flows  steadily  through  the  device  and  provides  the  downstream  gas  for 
entrainment.  The  secondary  manifold  uses  a  swirl-injected  gas  that  has  a  high  azimuthal  gas  velocity 
(relative  to  axial  velocity)  to  lengthen  its  axial  propagation.  This  swirled  gas  is  used  as  the  RMF-FRC 
formation  propellant.  In  this  design  the  slow  moving  gas  fills  the  small  end  of  the  cone,  while  the  fast 
moving  gas  fills  the  entire  cone.  The  RMF  ionizes  and  forms  an  FRC  in  the  small,  dense  conical  end.  It 
then  ejects  (or  is  pushed  with  secondary  acceleration  coils)  into  the  conical  acceleration  section, 
entraining  neutrals  and,  ideally,  lowering  total  effective  frozen  flow  losses. 
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Figure  12.  Photograph  of  the  PROTEAN  Figure  13.  Photograph  of  the  high  speed  SMBI. 
test  assembly. 


An  innovative  combined  circuit  has  been  designed  by  the  University  of  Washington  such  that  the 
same  accelerator  coil  can  be  used  for  both  a  bias  field,  as  well  as  pulse  high  field  acceleration.  This 
circuit  is  the  optimal  design  for  combining  an  external,  slow  bias  magnet  field  and  a  high  speed  pulsed 
accelerator  field,  all  in  the  same  magnetic  coil  and  with  minimum  charging  and  circuit  electronics.  This 
circuit  design  will  be  detailed  in  later  publications,  however,  the  magnetic  field  and  flux  profiles  given  here 
utilize  it.  Throughout  this  research  testing  was  initiated  with  vacuum  magnetic  field  testing.  The  circuit  was 
operated  first  with  only  bias  field  switches  and  was  found  to  be  well  within  design  tolerances.  Figure  14 
and  Figure  15  show  the  time  history  of  magnetic  field  and  magnetic  flux  for  coil  #1.  Figure  16  shows 
magnetic  flux  profiles  for  four  bias  field  configurations.  The  accelerator  field  coils  were  then  added  to  the 
circuit  with  magnetic  flux  measurements  for  all  the  first  accelerator  coil  group  shown  in  Figure  17.  Bias 
field  rise  times  is  set  so  that  the  RMF-formed  plasma  sees  a  constant  magnetic  field.  The  accelerator  rise 
time  is  chosen  such  that  magnetic  field  rise  times  corresponds  to  the  translation  time  of  the  magnetized 
plasma. 
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Figure  14.  Magnetic  field  history  for  the  bias 
field  discharge  in  coil  1.  Probe  is  located  4  cm 
from  the  base  of  the  discharge  insulator.  Shot 
#1073. 


Figure  15.  Magnetic  flux  history  for  the  bias  field 
discharge  in  coil  1 .  Probe  is  located  4  cm  from  the 
base  of  the  discharge  insulator.  Shot  #1073. 


Figure  16.  Composite  magnetic  field  profiles 
for  all  four  bias  field  discharges.  Axial  probe 
location  is  indicated.  Shot  #1 149. 


Figure  17.  Magnetic  flux  measurement  with  the 
primary  accelerator  coil.  Shown  is  5  V  bias  and 
100  V  acceleration  capacitor  banks. 


XENON  TESTING 


A  study  of  FRC  formation  using  a  xenon  propellant  and  various  propellant  injection  options  were 
investigated.  Using  only  the  annular  gas  injection  as  was  described  earlier,  a  robust  xenon  FRC  could  be 
formed  using  RMF  (Rotating  Magnetic  Field)  formation  at  1200  Volts.  This  is  a  total  of  2.9  Joules  and 
forms  a  well-coupled  FRC.  At  the  same  gas  pressure  and  puff  valve  settings,  neutral  gas  was  injected 
along  the  central  axis. 

In  addition  to  the  obvious  entrainment  effects  expected,  the  addition  of  neutral  gas  also  changes 
the  dynamic  behavior  of  FRC  formation.  The  residual  gas  on  axis  from  the  free-streaming  neutral  injection 
is  incorporated  into  the  RMF  FRC  and  acts  as  an  additional  source  propellant  for  formation.  As  a  result, 
for  the  same  injection  pressure  and  puff  timing,  the  resulting  FRC  appears  to  be  significantly  enhanced 
with  an  increased  ionization  rate  and  resulting  plasma  pressure.  Interestingly,  this  effect  was  not 
expected  based  off  of  larger  FRC  devices  tested  at  MSNW  which  optimized  with  an  annular  gas  injection. 
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In  order  to  match  the  plasma  conditions  and  FRC  source  conditions  as  closely  as  possible,  the 
injection  pressure  for  both  the  outer  annular  gas  injection  (to  be  called  RMF  propellant)  and  axial 
accelerated  neutral  propellant  (to  be  called  neutral  entrainment  propellant)  was  decreased.  The  results 
are  shown  in  Figure  18,  showing  for  all  conditions  well-formed  RMF  FRCs  with  rapid  ionization  and 
eventual  ejection  of  the  magnetized  plasma  are  formed.  When  gas  pressure  was  reduced  in  both 
propellant  injectors  to  15  psi,  a  similar  gas  condition  for  the  plasma  was  reached,  as  compared  with  no 
injected  gas.  While  identical  diamagnetic  signals  were  not  found,  this  condition  matched  appropriately  the 
plasma  formation  rate  of  ionization,  growth  within  the  discharge  chamber  region,  and  was  similar  in 
magnitude. 

With  a  similarly  matched  formation  condition,  the  downstream  effects  could  now  be  analyzed. 
Shown  in  Figure  19  is  the  approximate  plasma  density  as  determined  from  an  asymmetric  double  probe 
in  ion  saturation  in  supersonic  flowing  plasma.  As  is  typical,  the  absolute  error  of  these  probes  is  quite 
large,  though  it  is  believed  the  relative  error  between  operating  conditions  is  within  a  few  percent  for 
variability  for  velocity  and  peak  density  were  within  10%  for  these  measurements.  For  clarity,  error  bars 
will  not  be  shown,  as  they  are  not  critical  to  this  discussion. 

Downstream  results  are  quite  clear  and  shown  in  Figure  19.  When  there  is  no  neutral  injection 
and  RMF  FRC  is  formed  and  accelerates  through  a  constant  flux  expansion  into  the  outer  drift  chamber. 
The  ion  velocity  and  density  scale  as  expected  with  energy  and  the  RMF  propellant  injected.  With  neutral 
entrainment  injection,  a  similar  FRC  is  formed  and  injected  from  the  formation  section.  It  then,  however, 
quickly  decelerates  to  a  much  slower  speed  arriving  at  the  probe  well  over  100  microseconds  later. 

Radial  profiles  of  plasma  density  for  the  first  downstream  probe  are  shown  in  Figure  20.  Figure  20 
shows  a  highly  expanded  FRC  that  was  formed  in  a  4  cm  radius  cone  expanded  through  a  7  cm  radius 
cone  and  is  drifting  in  a  14  cm  cylindrical  section.  Plasma  densities  are  as  expected  for  this  magnitude  of 
expansion.  Interestingly,  the  expected  radial  profile  for  plasma  is  seen  with  a  lower  density  inner  core  and 
an  approximate  plasma  radius  of  8  cm.  While  this  lower  internal  density  is  expected  from  RMF-formed 
FRCs,  it  is  rarely  seen  downstream  due  to  the  relaxation,  inward  radial  ion  flow,  and  overall  geometry 
expansion  through  a  conical  magnetic  geometry.  This  annular  geometry  is  most  likely  due  to  the 
significant  annular  gas  injection  and  may  suggest  the  reasoning  for  the  increased  performance  with  a 
more  uniform  gas  distribution.  Additionally,  the  FRC  appears  to  be  a  coherent  magnetized  plasma  object 
moving  with  a  primary  bulk  velocity  (shown  by  arrival  time)  as  expected.  A  long  exposure  (250  ms) 
photograph  of  a  xenon  RMF  FRC  is  shown  in  Figure  21 . 


Figure  18.  RMF  Formation  with  and  without 
Neutral  Entrainment  propellant  injection.  Shot 
#’s  as  shown. 


Figure  19.  Downstream  double  probe  results  for 
the  key  RMF  formation  conditions. 
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Figure  20.  RMF  Radial  profile.  Shot  #1770-1842. 
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Figure  21.  Photograph  of  an  RMF-xenon 
discharge  at  approximately  2.8  Joules.  Annular 
pre-ionization  and  localized  ionization  region  is 
clear. 


NEUTRAL  ENTRAINMENT  WITH  XENON 

With  the  addition  of  the  downstream  on-axis  neutral  gas,  shown  in  Figure  19,  several  key 
differences  can  be  noted.  Firstly,  clear  deceleration  of  the  entire  plasma  body  is  seen.  In  addition,  axial 
compression  as  well  as  increased  radial  deceleration  and  radial  expansion  are  all  obvious.  These  are 
expected  effects  as  seen  in  the  previous,  referenced  modeling  efforts.  This  is  markedly  different  from  the 
fully  ionized  plasma  body.  A  radial  profile  of  the  resulting  plasma  configuration  is  shown  in  Figure  22.  The 
test  geometry  is  shown  again  in  Figure  23. 

Additionally,  the  peak  densities  of  the  decelerated  bodies  are  similar  to  the  initial  FRC  and  well 
within  the  probe  error  for  supersonic  plasma.  Effects  can  be  shown  quite  clearly  that  the  xenon  FRC 
appears  to  ionize  less  than  what  was  seen  in  the  neutral  neon  FRC  studies  performed  earlier.  Initial 
indications  are  that  this  is  due  to  the  lower  electron  temperature  in  the  small-scale  thrusters  compared  to 
the  larger  device  tested  previously  [12].  Another  difference  with  the  neon  results  are  that  a  central  plasma 
density  peak  approximately  2  cm  wide  is  evident,  which  is  lower  density,  but  high  velocity.  Rather  than 
toroidal  or  ellipsoidal,  the  overall  profile  is  wedge-shaped.  This  was  not  observed  in  either  the  modeling  or 
previous  neon  results.  Additionally,  the  peak  appears  to  have  drifted  off  axis  from  the  plasma  conditions. 
All  of  these  results  suggest  the  critical  dependence  on  injected  neutral  gas  and  that  this  device  would 
optimize  towards  an  even  higher  density  on  axis  that  previously  predicted  or  required  for  the  lighter 
gases.  Obvious  entrainment  and  deceleration  of  the  entire  bulk  plasma  is  seen. 
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Figure  22.  Radial  plasma  profile  for  an  RMF- 
formed  xenon  FRC  interacting  and 
decelerating  in  a  passive,  neutral  background 
gas.  Shot  #1770-1842. 


Figure  23.  Geometry  of  test  setup.  Shown  are 
Fast  Ion  Gauge  locations  1,2  and  Downstream 
TOF  Langmuir  probe  locations  1  and  2. 


Figure  25  shows  a  long  exposure  photo  of  an  RMF-formed  FRC  translating  into  a  background 
neutral  gas  region.  Some  increase  in  ionization  (bright  optical  emission  shows  excitation  during  the 
ionization  process)  during  the  downstream  region  underneath  the  first  acceleration  coil.  As  can  be  seen  in 
Figure  19  and  Figure  24  the  peak  arrival  of  a  plasma  body  (which  is  not  necessarily  and  FRC  at  this 
location)  is  much  later  than  any  simple,  passively  ejected  RMF  FRC.  Additionally,  compared  to  the  same 
diamagnetic  formation  signal,  the  peak  density  is  approximately  half.  It  is  also  worth  noting  that  there  is 
an  initial  high  velocity  peak  and  then  a  trailing,  slower  plasma  mass.  This  effect  will  be  discussed  in 
further  detail  in  the  later  acceleration  data. 


Figure  24.  Downstream  double  probe  results 
for  the  key  RMF  translation  conditions.  Shown 
are  two  RMF  formation  pressures,  downstream 
central  peak  arrival  with  a  neutral  gas,  and 
downstream  central  peak  arrive  with  the  addition 
of  acceleration. 


Figure  25.  Photograph  of  an  RMF-xenon 
discharge  with  a  downstream  neutral  gas.  Minor 
downstream  ionization  is  seen. 
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XENON  ENTRAINMENT  AND  ACCELERATION 


Acceleration  studies  were  performed  in  order  to  demonstrate  neutral  entrainment.  Ideally,  the  key 
metric  for  determining  the  success  of  the  neutral  entrainment  technique  was  the  ability  to  both  decelerate 
an  FRC  using  downstream  neutral  gas  as  well  as  accelerate  or  maintain  its  kinetic  energy  while 
increasing  the  available  propellant.  A  variety  of  coil  timing  tests  were  performed.  The  first  was  a  range  of 
empirical  timing  studies  for  the  first  acceleration  coil.  This  coil  is  situated  right  after  the  RMF  formation 
region  in  order  to  eject  the  FRC  at  optimal  temperature,  density  and  energy  without  having  to  rely  on  the 
passive  injection  methods  that  have  been  performed  in  the  past. 

This  study  varied  initial  acceleration  times  for  a  moderate  accelerating  field.  In  these  cases,  a  2.4 
kGauss  field  was  applied  as  shown  in  the  magnetic  field  results  in  Figure  26.  The  field  has  a  rise  time  of 
approximately  1 .5  microseconds  and  then  was  crowbarred  to  prevent  field  oscillation.  Empirically,  it  was 
found  that  the  optimal  injection  time  to  maximize  final  FRC  velocity,  as  well  as  final  diamagnetic  content 
(the  amount  of  excluded  field  from  the  FRC  which  thus  implies  the  magnetized  plasma  inventory), was  18 
+/-  2  microseconds  after  discharge  initiation.  This,  as  could  be  expected,  optimal  timing  occurred  at  the 
maximum  diamagnetic  and  plasma  conditions.  These  results  are  not  surprising  and  were  as  predicted 
and  discussed  earlier. 

However,  a  more  interesting  study  was  then  conducted  to  add  the  3  further  acceleration  coils.  A 
range  of  optimization  studies  were  performed  on  the  subsequent  coils.  For  the  results  shown  here,  a 
linear  separation  is  given,  though  MHD  results  show  that  for  an  accelerating  FRC  a  more  careful 
optimization  is  required.  For  so-called  constant  velocity  acceleration  in  which  mass  is  added  at  each 
stage  as  is  expected  in  the  steady  flowing  conditions,  a  linear  separation  is  likely  ideal.  In  this  linear 
operation,  each  subsequent  coil  is  triggered  fiber-optically  2  to  20  microseconds  apart  and  after  the 
primary  discharge.  Results  are  shown  in  Figure  27  and  several  key  scaling  can  be  found,  also  expected. 
First,  any  acceleration  appears  to  significantly  accelerate  the  FRC,  either  through  magnetic  gradient  or 
heat  addition.  If  the  timings  are  too  fast,  then  the  FRC  does  not  eject  from  the  device  before  being 
compressed.  Essentially,  instead  of  being  pushed,  it  is  crushed.  As  the  FRC  is  a  high-Beta  plasma 
configuration,  the  resulting  plasma  is  super-heated  and  compresses  radially,  increasing  radial  pressure.  It 
then  flows  along  axis  and  has  a  highly-centralized  radial  profile.  In  this  case  (5  microseconds),  the 
resulting  plasma  has  a  very  deformed  downstream  axial  profile,  appears  highly  centered  on  axis,  is  slow- 
moving  and  clearly  not  ideal.  As  the  coil  peak  field  timing  then  begins  to  follow  the  plasma,  the  FRC  found 
downstream  then  accelerates  and  the  plasma  seen  at  the  downstream  double  probe  arriving  sooner.  The 
timing  used  here  was  found  to  be  the  maximum  timing  tested,  trailing  the  FRC  significantly.  This  implies 
that  the  FRC  is  rather  insensitive  to  the  discharge  timing  as  long  as  it  is  downstream  of  the  peak  field.  It  is 
important  to  note  that  as  described  previously,  this  process  is  unlike  a  current-sheet  conical  acceleration 
as  the  plasma  diamagnetic  pressure,  and  thus  the  increase  in  applied  magnetic  field  gradient,  is  the 
primary  energy  transfer  mechanism,  rather  than  the  simple  cone  angle  of  the  acceleration  field  coils. 

Figure  28  shows  a  radial  profile  of  the  20  microsecond  case.  Figure  29  shows  a  250  ms  exposure 
photograph  of  an  RMF-formed  FRC  with  acceleration  through  a  background  gas.  No  major  increase  in 
ionization  is  seen.  The  radial  profile  has  several  key  elements  of  interest  at  this  stage.  First,  the  overall 
FRC  appears  to  have  expanded  further  radially.  This  was  somewhat  unexpected  as  in  the  MOQUI  MHD 
modeling,  the  FRC  was  always  shown  to  have  contracted  radially.  Additionally,  the  expanded  radial 
plasma  has  been  greatly  accelerated  (compared  to  the  non-accelerated  radial  profile).  This  expansion  will 
have  to  be  accounted  for  in  later  bias  field  design  to  force  radial  plasma  contraction  away  from  the 
thruster  walls.  Also,  the  central,  wedge-shaped  peaked  profile  has  been  reduced,  or  rather  the  large 
radial  extents  have  been  accelerated  more  than  the  center  region.  This  is  to  be  expected  if  there  is  less 
neutral  density  near  the  edges  and  unfortunately  makes  time  of  flight  velocity  measurements  quite 
problematic.  Perhaps  most  interesting  is  the  clear  transition  from  a  single,  peaked  profile  with  just  RMF  to 
more  ‘spread-out’  axial  profile.  There  is  a  clear  fast  moving  bulk  propellant  that  is  leaving  behind  a  tail  of 
slower  moving  plasma.  This  was  seen  somewhat  in  the  modeling  as  the  propellant  that  was  lost  from  the 
closed  field  region.  It  is  also  possible  this  is  an  artifact  of  the  pulsed  nature  of  the  neutral  injection  as 
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higher  density  neutral  gas  follows  the  FRC  out  of  the  puff  valve.  Follow-on  programs  will  seek  to  validate 
this  slow  moving  plasma  tail  with  a  steady  flowing  gas.  This  tail  will  have  ramifications  to  overall  thruster 
performance  and  gas  utilization,  if  realized  in  the  steady  flowing  case.  Rep-rate  and  overall  thruster  gas 
density  will  have  to  be  optimized  to  maximize  sweep-up. 


Figure  26.  Acceleration  magnetic  field  profile 
for  a  single  accelerated  FRC.  Total  energy  of 
0.75  Joules. 
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Figure  27.  Downstream  double  probe  results  for 
the  key  acceleration  conditions.  Shown  is  obvious 
acceleration,  over  compression,  and  no  acceleration 
conditions. 


Figure  28.  Radial  plasma  density  profile  for  an 
accelerated,  entrained  plasma.  Shot  #1770- 
1842. 


Figure  29.  Photograph  of  an  RMF-xenon 
discharge  with  a  downstream  neutral  gas  and  2.3 
Joules  of  acceleration. 
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CONCLUSION 


A  new  innovation,  the  PROTEAN,  adds  a  secondary  Neutral  Entrainment  (NE)  stage  that,  if 
succesfull,  can  reduce  plasma  formation  losses  in  electromagnetic  and  electrostatic  electric  propulsion 
thrusters.  This  NE  stage  injects  neutral  propellant  into  the  plasma  path  after  ejection  from  an 
electromagnetic  Field  Reversed  Configuration  (FRC)  thruster,  such  as  the  1  kW  Electromagnetic 
Plasmoid  Thruster  (EMPT).  The  FRC  ingests  these  neutral  particles,  accelerating  them  up  to  high  velocity 
through  resonant  charge  exchange  collisions.  The  NE  stage  then  efficiently  adds  further  kinetic  energy  to 
the  engorged  plasmoid  through  Peristaltic  Dynamic  Acceleration,  in  which  sequenced  theta-pinch  coils 
are  actuated  to  provide  a  large  magnetic  pressure  gradient.  Presented  are  results  from  a  single-pulse 
experiment  which  show  that  Neutral  Entrainment  works  as  theroetically  predicted  for  a  xenon  propellant. 
Xenon  RMF-FRCs  were  formed,  injected  into  a  neutral  gas,  and  then  accelerated  through  that  neutral 
gas.  A  wide  range  of  timing  and  gas  optimization  tests  were  done.  They  showed  clearly  that  the  xenon 
FRCs  can  be  decelerated  significantly  by  the  centrally-located  neutral  gas  and  then  with  the  application  of 
four  pulsed  magnetic  field  coils,  can  be  accelerated  back  up  to  high  velocity.  This  appeared  to  happen 
without  wall  interaction,  with  minor  additional  ionization,  and  no  apparent  increase  in  frozen  flow  losses. 
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msnw*  Electric  Propulsion  on  Light  Gases 


Real  plasma-based  propulsion  systems  take  5-1  OX  minimum 
ionization  energy  to  actually  form  and  eject  a  plasma 


* 


Electrostatic  and  Electromagnetic  systems  pay  a  large  efficiency 
penalty  to  operate  on  lightweight  gases  (also  T/P  for  heavy  gases) 


Specific  Impulse  [si 

Theoretical  maximum  thruster  efficiency  with  a  40  eV/ion 
frozen  flow  (plasma  formation)  loss 


How  can  we  accelerate  lightweight,  chemically  reactive  gases,  but  without  the  large 

ionization  costs? 
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&>i  msnw*  Neutral  Entrainment 

ISSUE:  Plasma  formation  is  the  primary  loss  mechanism  for  ALL 
electrostatic  or  electromagnetic  propulsion  systems.  It  prevents 
operation  at  lower  specific  impulse,  lowers  efficiency  at  all  exit  velocities, 
and  requires  high  mass  propellants. 

Solution:  Entrain  neutrals  in  an  acceleration  field  after  formation. 


Benefit:  Specific  impulse  can  be  tailored  to 
the  mission,  while  the  thruster  operated 
at  its  maximum  efficiency,  even  on  light 
propellants. 

•  Basic  Idea:  An  FRC  will  ‘injest’  large 
quantities  of  neutral  gas  through  charge 
exchange  collisions,  not  ionization.  If 
you  accelerate  an  FRC  while  providing 
upstream  neutral  gas,  Isp  can  be 
specified  and  mass/thrust  added  with 
very  high  efficiency  [1]. 


[1]  Matsuzawa,  Y.,  et.  al,  “Effects  of  background  neutral  particles  on  a  field -reversed 
configuration  plasma  in  the  translation  process”.  Phys.  Plasmas  15,  (2008). 


Dynamic  formation  and  acceleration  of  an  FRC,  shown  with  4  field  coils 
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Key  Concept 


•A  MSNW> 


Entrainment  of  Neutrals  in  a  Magnetized  Plasma 


■=> 


High  energy 
ions  enter 
neutral  gas 


Plasmoid  ions 
charge  exchange 
with  neutral  gas 


New  magnetized  ion 
confined  by  magnetic 
field 


New  magnetized  ion 
confined  by  magnetic 
field 


•  Charge  exchange  collisions  dominate 

•  Ionization  must  be  minimal 

•  Gas  utilization  is  key 


Effectively  accelerate  neutrals  rather  than  ions 
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msnw>  Field  Reversed  Configuration  (FRC)  /i 

Magnetic  Field  lines  and  Pressure  Contours  V 

fs'j 


z 

Flux  Conservation 

External  measurements  of  B  yield 
FRC  separatrix  radius  rs(z),  FRC  length  Ls 
=^>  volume,  position,  velocity 

Radial  Pressure  Balance 

Simple  cross-tube  interferometric 
measurement  with  rs  from  yields  <n>  and  T 

Axial  Pressure  Balance 

With  above  obtain  plasma  energy, 
Inventory,  confinement  times 

FRC  equilibrium  constraints  and  the  diagnostic  measurements  that 
together  with  the  equilibrium  relations  that  are  employed  to  determine 
the  basic  parameters  of  the  FRC  equilibrium 
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Equilibrium  ^ 
Relations: 
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MSNW> 


2D  Resistive  Calculation  of  Acceleration 


Standard  FRC  Thruster  Operation 


FRC  thruster  with  sequenced  coil  array 


t  =  8.01  |isec 


“Afterburner”  following  current  FRC  thruster  accelerates  FRC  through  peristaltic 

electromagnetic  acceleration 
Increase  total  dM/dx  through  multiple  field  coils 
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Ion  -  Neutral  Interaction 

Employing  the  SEL  Code 


-  SEL-HiFi  Visco-Resistive  MHD 


-  Neutral  and  Ion  Fluids  added 

-  Flux-source  form  for  inputs 

-  Completely  implicit  with  an  adaptive  grid  with  high-order  spectral  elements 

-  Neutral  fluid  with  resonant  charge  exchange,  electron-impact  ionization,  and  radiative 
recombination  reactions 


Pn»  k8/m3 
2.500e-05 
1.875e-05 
1.250e-05 
6.250e-06 
0.000 


t=9  jis 


t=27  \is 


t=36  |is 


Neutral  density  psuedocolor  plots  for  baseline  ELF  simulation  with  a  Gaussian  radial  neutral  gas  profile.  Black  and  white  arrows  indicate  the  direction  and 
magnitude  of  plasma  and  neutral  momentum,  respectively.  Clear  neutral  entrainment  is  seen. 


Meier,  E.T.,  et  al.  "Development  and  validation  of  a  two-fluid  plasma-neutral  model",  Innovative  Confinement  Concepts  (2011). 
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axial  momentum,  kg  m/s 


Pressure  psuedocolor  plots  for  baseline  ELF  simulation.  Clear  neutral  entrainment,  pulse  sharpening,  and 
increase  in  external  magnetic  field  pressure  are  seen. 


xIO?  &d0'3 


Temporal  profiles  of  momentum-related  quantities  for  ELF  baseline  simulation 
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Experimental  Program 


1 .  Neutral  Entrainment  demonstrated  as  part  of  the  AFOSR  MACEEP  pro^tci'lfi 

2.  Phase  I  of  PROTEAN  to  demonstrate  neutral  entrainment  in  a  thruster-like 
environment 

3.  Demonstrate  entrainment,  deceleration,  and  magnetic  acceleration 


Time 


z  (cm)  PROTEAN  Phase  I.  Shown  is  a  Xenon  FRC. 

MOQUI  modeling  of  a  Conical  acceleration  section  with  an  RMF  FRC 

SCIENTIFIC  CHALLENGES: 

•  Predict  and  experimentally  validate  the  velocity  and  density  limitations  of  neutral  entrainment  in  a  thruster  application 

*  Quantify  effects  of  charge  exchange  and  ionization  with  an  easily  ionized  and  thruster  relevant  gas  (Xenon) 
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Neutral  Injector 


jt 


>  Neutrals  injected  on  axis  and  in  limited  axial  location  to  isolate  effects 

>  Supersonic  Molecular  Beam  Injector  designed  and  implemented 

>  Clear  beam  formed  with  +/-  3  degree  ejection  angle 

>  Coaxial  neutral  beam  injector  and  tangential  swirl  pre-ionization  injector 


-5 


Time[s] 


x  10 


Coaxial  SMBI  and  pre-ionization  electrode 


Neutral  Density  profile  in  mT  for  the  discharge  of  interest  showing 
initial  neutral  beam  and  follow-on  FRC  propellant  flows  at  68  cm. 
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{<♦3 msnw>  PROTEAN  Circuit  and  Accelerator 

•  Combined  bias  and  accelerator  circuits 

•  Capable  of  simple  energy  recovery 

•  5  total  bias  field  an  accelerator  coils 


soce 


•  50  -  200  G,  1 0  ms  bias  with  500  -  3000  G,  3  jis  accelerators 


Magnetic  flux  measurement  with  the  primary  accelerator  coil. 

Shown  is  5  V  bias  and  100  V  acceleration  capacitor  banks. 


Composite  magnetic  field  profiles  for  all  four  bias 
field  discharges.  Axial  probe  location  is  indicated 
Shot  #1149. 
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Xenon  FRC  Formation 


£3  msnw*  Xenon  FRC  Formation 

2.8  J,  300  kHz  RMF  formation 

FRC  formed  and  passively  translated  downstream 

Expected  downstream  velocities  and  arrival  times  measured  with  TOF  double 
Langmuir  probes 

Gas  only  injected  with  outer  formation  injector  (swirl) 


e"ce 


Photograph  of  an  RMF-xenon  discharge  at 
approximately  2.8  Joules.  Annular  pre-ionization 
and  localized  ionization  region  is  clear. 
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Internal  neutral  beam  added 


Xenon  Entrainment 


e*ce 


Downstream  profiles  matched  (as  much  as  possible)  between  only  formation  gas 
injection  (outer  swirl)  and  inner  neutral  beam  injection 

Clear  deceleration  and  apparent  entrainment 

No  increase  in  downstream  density  or  optical  radiation  (ionization) 

Radial  and  axial  expansion  and  ‘smearing’  of  FRC  profile 


Geometry  of  test  setup.  Shown  are  Fast  Ion  Gauge  locations 
1,2  and  Downstream  TOF  Langmuir  probe  locations  1  and  2 


Downstream  double  probe  results  for  the  key  RMF 
translation  conditions.  Shown  are  two  RMF 

formation  pressures,  downstream  central  peak 
arrival  with  a  neutral  gas,  and  downstream  central 
peak  arrive  with  the  addition  of  acceleration. 
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^ MSNW*  Xenon  Entrainment  and  Acceleration  /Trr\ 

•  Acceleration  fields  added  to  decelerated  profile  -  2.4  kG,  2.3  Joules  v  yS 

•  Single  initial  acceleration  in  formation  region  optimized  18  microseconds  after  f>SICen0 
discharge  iniation  (slightly  past  peak  diamagnetic  signal). 

•  3  secondary  acceleration  coils  added  with  uniform  spacing  (5-20  microseconds) 

•  FRC  could  be  accelerated  back  up  to  initial,  non-entrained  velocity 

•  Optimal  timing  was  surprisingly  late,  i.e.  peak  field  always  trails  peak  diamagnetic 
signal  (both  geometrically  and  temporally) 

•  Apparent  increase  in  slow  moving  tail  plasma,  transient  gas  injection  is  likely  culprit 


Radial  plasma  density  profile  for  an  accelerated,  entrained 
plasma.  Shot  #1770-1842. 


Downstream  double  probe  results  for  the  key 
acceleration  conditions.  Shown  is  obvious 
acceleration,  over  compression,  and  no  acceleration 

conditions. 
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to  msnw*  Conclusions  Ai 

ff  mi 

•  Neutral  Entrainment  has  the  capability  to  reduce  effective  plasma  frozen  flow 

losses  in  electromagnetic  propulsion  systems  Ic 

•  The  PROTEAN  test  bed  was  built  and  demonstrated  Xenon  and  Argon  FRCs 

•  Combined  pulsed  magnetic  bias  fields  with  peristaltic  acceleration  demonstrated 

•  Xenon  operated  with  low-pressure  (temperature)  FRCs  showed  little  ionization 

•  Transient  gas  injection  lead  to  radial  and  axial  expansion  and  gas  under-utilization 

•  Xenon  RMF-FRCs  were  formed,  injected  into  a  neutral  gas,  and  then  accelerated. 

•  Xenon  FRCs  can  be  decelerated  by  the  centrally-located  neutral  gas  and  with  the 
four  pulsed  magnetic  field  coils,  can  be  accelerated  back  up  to  high  velocity. 


Downstream  double  probe  results  for  the  key  acceleration 
conditions.  Shown  is  obvious  acceleration,  over 
compression,  and  no  acceleration  conditions. 
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Photograph  of  an  RMF-xenon  discharge  with  a  downstream 
neutral  gas  and  2.3  Joules  of  acceleration. 


Distribution  A:  Approved  for  public  release;  distribution  is  unlimited. 


BACKUP -Neutral 
Entrainment 
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MSNW> 

Field  Reversed  Configuration  (FRC) 

Magnetic  Field  lines  and  Pressure  Contours 


Key 

Equilibrium  < 
Relations: 


P0=n0kT=|^ 

2no 


Flux  Conservation 

External  measurements  of  B  yield 
FRC  separatrix  radius  rs(z),  FRC  length  Ls 
=^>  volume,  position,  velocity 

Radial  Pressure  Balance 

Simple  cross-tube  interferometric 
measurement  with  rs  from  yields  <n>  and  T 


Axial  Pressure  Balance 

With  above  obtain  plasma  energy, 
Inventory,  confinement  times 


FRC  equilibrium  constraints  and  the  diagnostic  measurements  that 
together  with  the  equilibrium  relations  that  are  employed  to  determine 
the  basic  parameters  qf.the.FRC,  equilibrium  ,  ,  „ 
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MSNW> 

Increase  in  FRC  Isp  with  Follow-on 
Magnetic  Accelerator 

80 

60 

Vfrc 

(km/s)  40 
20 
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0  4  8  12  16  20 

Time  (psec) 

Kinetic  energy  gain  observed  with  accelerator  would  be  employed  to 
maintain  FRC  Isp  during  entrainment  i.e. 

F  =  mVFRC  rather  than  F  =  mVFRC 
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2D  Numerical  Calculation  of  Ion  -  Neutral  Interaction  * 
Employing  the  SEL  Code  -  Initial  conditions  PSI 

D  Ion  Density 
(xIO19  rrr3) 

3.000 
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—  0.7500 


0.000 


D  Neutral  Density 
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2D  Resistive  Calculation  for  Entrainmeritr7 

1‘SfCentef 

Experiments  with  ELF  thruster 


Standard  ELF  Thruster  Operation 


ELF  thruster  with  sequenced  coil  array 

t  =  8.01  |isec 


“Afterburner”  following  current  ELF  thruster  would  be  used  to  increase  thrust  at 
constant  Isp  at  very  high  efficiency  with  neutral  entrainment 
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SEL  Implementation  of  Visco-Resistive  MHD 
with  a  Dynamic  Neutral  Fluid  Background 


Density 

Equations 


dp, 


dt 

dp„ 

dt 


+v.(Avi)  =  H/i(r,i)„-r,,f) 


+v*(A,v„)=mj(rref-rj(„) 


5 


T-( A V, )  +  V*(p,  v(  v,  +  pi)  =  j X B  - r recm,\,  +  r h„m,vn  +  r„  ( m, v„  - mf v, ) 

Momentum 

Equations  JL(  p,v  „ )+ V(p„v„v„ + Pni) = r„fw,v,  -r  +r„  (m,v,  -  m,\„ ) 


1  dp 
y-\  dt 

Energy 

Equations 


dt 


+  V< 
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j-i 


pyi 


=  v,*P  +  (  r„  +rion)(KE„  +  KEi-miv,yn)  +  rion(KEn-<fiion) 


1 


/-I  dt 


+  V. 


r 


r-i 


=v„*vp„ 


+  r,,c  (KE,  +  KE.  -  $ )  +  r„  ( AE,  +  KE„  -  miv„ -v, ) 


•  Single-fluid  plasma  model  which  allows  for  Spitzer-Chodura  resistivity, 
Braginskii  thermal  conduction,  symmetric  viscosity,  and  dissipation 

•  For  readability,  dissipative  terms  are  excluded. 
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Atomic  Physics  Rates 


The  equations  being  used  for  atomic  physics  rates  are 

T  =2*54.82*10“20T  -2057  * 
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where  the  ion-neutral  relative  velocity,  vm,  is  computed  from 


( 


= 


A 


1/2 


V  mi 


kT. 


J 


m 

s 


•  It  is  assumed  that  the  radiation  after  recombination,  hv,  is  immediately  released 
and  is  equal  to  cpion.  It  is  also  assumed  that  rrij  =  mn 

•  The  equation  for  charge  exchange  is  based  on  a  curve  fit  of  data  in  the  1966 
NASA  Technical  Note  by  Nichols  and  Witteborn 
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>  The  KEY  idea:  A  price  in  efficiency  MUST  be  paid  to  operate  in  a  lower  mass  gas. 
Fundamental  physics. 

>  There  are  desirable  mission  reasons  to  be  able  to  use  air  gases,  i.e.  no  more  rocket  equation 

Assumptions  in  Model  and  Previous  Graphs 

•  Frozen  flow  losses  are  conservative: 

-  40  eV/ion  plasma  formation,  radiation,  ionization,  wall  losses,  etc.  (worse  than  measured) 

1 0  eV  total  residual  plasma  temperature  (Te  +  Ti).  (worse  than  measured) 

•  System-level  losses: 

-  95%  PPU  efficiency  (better  than  measured,  but  demonstrated  elsewhere) 

-  95%  Pulsed  inductive  energy  addition  with  IGBT  switching  (not  demonstrated,  but  modeled) 

•  Neutral  entrainment  CE  efficiency  : 

1 00  %  Resonant  Charge  Exchange  (small  changes  would  be  expected  for  non-resonant  charge  exchange  with  02  to  be 
explored  via  model) 

-  N  is  number  of  charge  exchanges  per  initial  FRC  ion 

•  Current  Thruster  Data 

Uncertainties  are  both  shot-to-shot  and  absolute  measurement  uncertainty 

-  Thruster  efficiency  does  not  include  PPU  losses  and  uses  measured  downstream  velocity  and  impulse 
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